Cavity-enhanced detection of magnetic orders in lattice spin models 
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We develop a general scheme for detecting spin correlations inside a two-component lattice gas of 
bosonic atoms, stimulated by the recent theoretical and experimental advances on analogous systems 
for a single component quantum gas. Within a linearized theory for the transmission spectra of the 
cavity mode field, different magnetic phases of a two-component (spin 1/2) lattice bosons become 
clearly distinguishable. In the Mott-insulating (MI) state with unit filling for the two-component 
lattice bosons, three different phases: antiferromagnetic, ferromagnetic, and the XY phases are 
found to be associated with drastically different cavity photon numbers. Our suggested study can 
be straightforwardly implemented with current cold atom experiments. 
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I. INTRODUCTION 

Atomic quantum gases trapped in optical standing 
waves have become ideal systems for implementing lat- 
tice spin models after the pioneering theoretical proposal 
[l[ and the experimental observation Q of the super- 
fluid (SF) to Mott insulator (MI) transition in the Bose- 
Hubbard model. When atoms of two-species or two- 
components are loaded into an optical lattice, a vari- 
ety of more general effective spin models can be con- 
structed [H, U [H, including the well-known anisotropic 
Heisenberg XXZ model. The development of noise spec- 
troscopy @, 0, S, ©, M, El has provided an astounding 
breakthrough that overcomes several significant hurdles 
in detecting quantum correlations, or in measuring the 
second order spin moments for various magnetic phases 
of lattice models. 

Cold atoms are usually probed with time of flight 
methods, which measures the atomic density or matter- 
wave interference patterns upon being released from 
traps and often after significant expansions. The near 
resonant imaging light generally destroys the atomic 
state. Several quantum limited detection schemes have 
since been suggested, capable of quantum non-demolition 
detections of strongly correlated states in atomic lattice 
models [HI, [l3j |. A very interesting approach relies on 
the enhanced detection sensitivity provided by an optical 
cavity, as was first proposed by Mekhov et, al. 0, [ll| . 
The transmission spectra, calculated to the first order, 
or within the linear response theory of the amplitude for 
the probe field, assumes the initial state of atoms to re- 
main unchanged when expectation values are taken and 
carries unambiguous signatures of magnetic orders in an 
atomic Bose-Hubbard model. 

Several experimental groups have recently succeeded in 
the difficult first step of coupling atomic condensates into 
high Q optical cavities [m[l7f. highlighting the prospects 
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FIG. 1: (Color online) Schematic illustration of the proposed 
experimental setup and the level diagram for a bosonic atom 
with two states resonantly coupled to the two cavities. 



for creating and detecting exotic quantum phases of lat- 
tice spin models [18|. A promising new direction worthy 
of theoretical investigation concerns the study of atomic 
lattice spin models coupled with optical cavities, gener- 
alizing the single component study [13, [l5j]. Nonlocal 
quantum spin correlations of the various magnetic orders 
could analogously be reflected through the photon num- 
bers and statistics. 

This paper describes a scheme for detecting spin corre- 
lations in a two-species or two-component bosonic atom 
lattice [1, 0, 0, H3] ■ Our study shows that atomic spin 
correlations are faithfully mapped onto the transmission 
spectra of the cavity probe field, making them easily di- 
agnosed through cavity QED based techniques. 
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II. MODEL 

Our model is based on the scattering of two Raman 
matched incident laser beams from a lattice of effective 
spin 1/2 bosonic atoms [2l], [52]. Similar to the original 
model [l5f for single component bosons, we consider TV 
atoms with two internal states identically trapped in an 
optical lattice with M sites formed by far off-resonant 
standing-wave laser beams. As schematically illustrated 
in Fig. [T] K < M lattice sites are located within the over- 
lapped region of the two fundamental modes of the cavi- 
ties. We consider two non-degenerate hyperfine states |1) 
and 1 2), the two stable ground states that are coupled to 
a common excited state |3) with a blue common detuning 
A and no differential detuning, forming a Raman coupled 
A-type atom model. The resonant cavity modes are de- 
noted by matching labels with frequencies uii (I = 1,2). 
For large detu ning A, we adiabatically eliminate the ex- 
cited state 1 3} j24l |25|| and end up with two-state atoms 
effectively coupled in the overlapped region of two opti- 
cal cavities. For a single atom, the effective coupling is 
described by f2a{a2&{&2 + h.c. and the ac Stark shift be- 
comes 8ia,aib ,b a —i with Si = <7?/A and = gi^/A. 
The peak value for the dipole coupling with their respec- 
tive cavity mode is denoted by g a for transition between 
\a) «-» |3). 6<j = i.2 (0^=1,2) denotes the corresponding an- 
nihilation operator for the atom (cavity mode photon). 

Following the notations of Ref. [HI , the Hamiltonian 
for effective spin 1/2 bosons in a lattice coupled to two 
optical cavities takes the form Hb + Hj, with 

Hj = fuviajai - ihrj (aie tbJlpt - h.c.) 
1=1.2 

K K 

+ Mi Wlf "-il«l a l + ^2 ^2 I" 2 ! 2 n i2 a 2 a >2 

t=l i=l 
K 

+ hnJ2 (Aialazb^ba + h.c?j , (1) 
1=1 

where n J(T = &j .&i c7 gives the number of atoms in 
state \a) at site i and 01,2 (r) is the mode function 
of the cavity with wave- vector ki.2- The coefficients 
Ai{9\,6 2 ) = u*(ri) u 2(i"i) due to emission/absorption or 
absorption/emission cycle are responsible for the geomet- 
ric dependence of the effective coupling [151 ] . 

With atoms assumed to occupy only the lowest Bloch 
band, our model generalizes the familiar Bose-Hubbard 
for two-components: Hb as in Eq. (1) of Ref. [l!| for two 
species. Following the work of [la ], we perform a linear 
calculation to the first order in cavity probe field, thus we 
leave out the dynamics of how various quantum phases 
of the atomic lattice are realized or dynamically created 
through the tuning of lattice parameters. This further 
justifies the neglect of atomic tunneling as well as the 
on-site intra- and inter-component interactions. In addi- 
tion to the coupling of each atomic component with its 
corresponding cavity mode, Raman matched two-photon 
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TABLE I: Cavity 1 photon number for the four quantum 
phases of the two-component Bose-Hubbard model at the 
diffraction maxima (minima) with 6\ = (6\ = tt/2) and 
02 = 0. For the XY phase 6 A = 6b = tt/3. 

processes can transfer atoms between the two effective 
spin states, unless the atoms are prepared in the so-called 
dark state |dark) ~ (02)32 |1) — (o-i)gi |2) corresponding 
to Coherent Population Trapping (CPT) [Hj|. We also 
assumed large detuning between cavity and atoms, to 
keep the actual excitations low, or negligible; thus any 
Raman type population transfers only affect the initial 
state to higher orders than the linear response theory 
calculation we provide. The second term in Eq. (fTJ) de- 
scribes the coherent pumping of cavity 1 at frequency uj\ p 
with amplitude 77. 

III. SEMICLASSICAL THEORY 

We first consider the simplest case with no external 
pumping on cavity 1, i.e., r\ = 0, and assume cavity mode 
ct2 to be a classical field, or a c-number amplitude as in 
Ref. [l5l ]. In the frame rotating with frequency LO2, di 
evolves in time according to the Heisenberg equation 

K 

hi = — i(Ai2 + 61 y~] |ui| 2 n i i)ai - iflDa 2 - «ffli, (2) 

i 

where A12 = to\ — ^2 and k denotes the cavity decay 
rate and is put in by hand. Its corresponding Langevin 
noise is neglected. We have defined the analogous opera- 
tor D = Yli=i AiS~ , in terms of the effective lattice spin 

operators = &i2 and 5 2 + = (S^)\ which obey the 
standard commutation relation at the same site and com- 
mute with each other on different sites. Neglecting the 
presumably much smaller cavity field induced ac Stark 
shift in comparison to A 12 or k ai and the photon 
number is easily obtained as 

ai = CD, a\ ai = \C\ 2 D^D, (3) 

where C = — iVLa^/ii A12 + k). The photon number 
(a\ai), clearly provides information about the spin cor- 
relation in the two-component bose lattice through the 
moments associated with the same site (S^ Sf) and be- 
tween the different sites (S^S~). The angular depen- 
dence can become totally different due to the geometric 
coefficients Aj(#x,#2)- Within the linear response, the 
above averages are expectation values with respect to 
whatever initially prescribed atomic ground state. 
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The quantum phases for a two-component lattice 
bosons at commensurate fillings have attracted signifi- 
cant attention [l9|,[2(|. The phase diagram consists of (1), 
2MI where both boson components are in the MI phase; 
(2), SF+MI where one is SF and the other is MI; and (3), 
2SF where both components are SF. Deep inside the MI 
phase the ground state of the system may be character- 
ized by filling the lattice site with even or odd numbers 
of atoms [2(| • In addition to the usual even filling phase 
with tlx = n 2 , a particularly interesting phase arises when 
the total filling factor is odd, especially at unit filling, 
i.e., for ri\ + n 2 = 1. This exotic phase has been ex- 
tensively studied 0, BL [2 H3] by adopting a trial wave 
function \^mi) = ElieAjeB \^A)i\ipB)j, which is of a 
form composed of two sublattices A and B with \iPa,b) = 
cos(^,B/2)|l,0) + e i ^. B sin(^ A ,B/2)|0,l)- \ni,n 2 ji de- 
notes the state with n\ (n 2 ) number of component- 1 
(-2) atoms at site i and 9s and </>s are variational pa- 
rameters. Three types of spin exchange interactions 
are identified: (I), anti-ferromagnetic phase (AF) with 
9a = 0(7r) and 9b — tt(0); (II), ferromagnetic phase 
(FM) with 9 A = 9 B = 0; and (III), XY phase with 
9a = 9b 7^ 0. The 2SF phase, whose quantum state 

is *sf ~ (Ei&li^tEi^lO) with ^1,2 the total 
number of component-1 (-2) atoms [2(| . will serve as a 
reference for presenting our results. 

The scattered photons are explicitly tabulated in Ta- 
ble [U For a ID optical lattice of a spatial period d = A/2 
and with atoms trapped at sites centered at Xj = jd, the 
mode functions are ui^i^j) — exp^jjki^ldsinf?!^) for 
a traveling wave and/or u\ t 2 (jj) = cos(r/|ki. 2 |<i sin ^1,2) 
for a standing wave form. Atoms in the FM phase do 
not scatter because the two coupling paths to the ex- 
cited state 1 3} destructively cancels as in the dark state. 
For the notation we use, the FM state corresponds to 
all atoms staying in state |1), then a semi-classical light 
amplitude (02) g 2 clearly will not be able to cause any 
scattering. While the initial atomic states of the AF 
and XY phases under the single excitation of a semi- 
classical light are not any more dark states, they will 
scatter. These features thus completely characterize the 
many-body spin correlations of the quantum phases for 
the two-component Bose-Hubbard model. 

To map quantum fluctuations of lattice spins faith- 
fully onto the probe cavity photon statistics, we define 
a noise function i?(6q,6» 2 ) = (D^D) - {D^)(D), whose 
angular distribution is compared in Fig. [2] for all four 
quantum phases. The structure in the angular distribu- 
tion comes from the summation of the geometric coeffi- 
cients from different sites, reflecting both the on-site and 
off-site lattice spin correlations. In the SF phase with 
rii = n 2 = 1/2, the respective noise functions are com- 
pletely different for the two choices of cavity modes. For 
the traveling wave, the noise function is zero for the FM 
phase, but takes nonzero values and is isotropic for the 
XY and the AF phases. The angular dependence for the 
standing wave mode case is richer than that for the trav- 
eling waves. The structures in the angle dependence can 
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FIG. 2: (Color online) The angular distribution of 11(61,82) 
for the four quantum phases evaluated for different choices 
of cavity mode functions: the left (right) panels are for two 
traveling (standing) waves and for 82 — (82 = O.I71-). We 
have assumed N — M — 2K — 40 and in the SF phase 
m = n 2 = 1/2. For the XY phase 8 A = 8 B = tt/3. 



be attributed to dependence on the summation of the 
geometric coefficients, and physically due to both on-site 
and off-site lattice spin correlations. 



IV. QUANTIZED MODEL 

We next consider the more general case with coherent 
pumping for cavity 1 at frequency uj\ p [l5j . The dissi- 
pations for both cavities are assumed the same with the 
associated Langevin noise terms neglected in the Heisen- 
berg operator equations. Within a linearized calculation, 
we decorrelate the atomic and field operators and replace 
in the Heisenberg equations for a\.i the atomic opera- 
tors by their respective expectation values, which leads 
to (at )(a/) = |(a;}| 2 . To simplify our result, we further 
assume |ui.2(ri)| 2 = 1, which occurs for the diffraction 
maxima with A; = 1 at 9\ — or the minima with 
Ai = (—I) 1 at 9\ — ir/2 when the ID lattice is lined 
up at 02 = 0. The cavity photons are found to be 

(a\)( ai ) = r,V + (i)/B 1 (a\)(a 2 ) = i^otajB, (4) 

where B = n 4 + K 2 ((f + Cl + 2a* a) + (C1C2 - a* a) 2 , 



a = n Ef MSrh and Cl = A ip + 5 l Ef (n«). The de- 
tuning Ai p = ui — lu± p are assumed the same for I = 1 , 2 
because A12 <S u>i,2- If the cavity coupling is assumed 
identical, we end up with £1,2 = Cl = S. Equation (U|) 
shows that probe photon numbers depend on the aver- 
age values of on-site atom numbers (rii a ) and the lattice 
spin operators (5- ). A crucial term for spin correla- 
tion a* a appears in the expression for (a\)(a2)- At the 
diffraction minima or maxima a — so that no pho- 
ton will be detected from cavity 2 except for the XY 
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phase. This then allows for simplified expressions of the 
scattered photon numbers (a\)(ai) from the AF and FM 
phases into (a|)(a 1 ) = t] 2 /(k 2 + ( 2 ), which only depends 
on the detuning Ai p and atom numbers for component- 1 
in the overlapped if-sites = JY ( n n)- When a^O, 
however, (a\)(ai) for the XY phase at the diffraction 
maxima depends on two parameters Ci an d C2 includ- 
ing the detunings A lp , A 2p , and the number of atoms 
for both components in the overlapped region of if -sites. 
Measuring photon numbers (a\)(ai) thus gives sufficient 
information to distinguish magnetic orders or quantum 
phases of the two-component Bose-Hubbard model. 

An especially interesting property concerns the depen- 
dence of the probe photon numbers on the detuning Ai p , 
as is illustrated in Fig. for the four quantum phases. 
For the FM and AF phases, we find Lorentzians with 
width k and shifted by SN^ as in the classical result of a 
single component Bose-Hubbard model [l5j]. In contrast, 
for the SF phase the photon number distribution is an 
envelope of a comb for a good cavity (k = 0.1(5) while a 
smooth broadened contour for a bad cavity (k = S). In 
the SF case, individual atoms are completely delocalized 
over all sites causing significant number fluctuations over 
each site within the K-site region. The corresponding 
quantum state is a superposition of Fock states contain- 
ing all possible distributions of atoms for component - 
1 at if sites, which gives rise to scattering terms from all 
possible atomic distributions. For the XY phase, the dou- 
ble peaked feature provides evidence for different popu- 
lation of atoms in the two internal states, with the rel- 
ative heights of the two peaks being controlled by the 
variational parameters &a,b- This structure in the easy 
plane XY phase is essentially identified with the so-called 
superfluid counterflow (SCF) phase, which can be quali- 
tatively understood as a paired superfluid vacuum (PSF) 
phase, a strongly correlated superfluid ground state al- 
ready predicted from numerical simulations These 
distinct features of the transmission spectrum we discuss 
for the various quantum phases form the basis for easily 
detecting and differentiating the corresponding magnetic 
orders in the two-component Bose-Hubbard model. 

Like the original cavity scheme of Mekhov et. al. 
[3, UH, the scheme we propose, is constructed to detect 
high order moments. The different phases (in the sense 
of quantum states of matter) of a two-component lattice 
bose gas are resolved from the statistics of scattered pho- 
tons or pseudo-spins. In this sense, it is analogous to the 
so-called noise spectroscopy of quantum gases @, [28| , al- 
beit somewhat superior due to the enhanced collection 
efficiency aided by cavities. The Ramsey spectroscopy 
[27| . as proposed by Kuklov, measures the first order 



moments of atomic pseudo-spins. The SCF state or the 
paired condensation phase is a special case, where the 
order parameters are simply field operators themselves. 
Thus their presence can be probed by the Ramsey spec- 
troscopy measurement of the relative phase (in the sense 
of amplitude and phase) . 
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FIG. 3: (Color online) Cavity 1 photon numbers as a function 
of cavity-probe detuning for the four quantum phases: AF 
(red dashed dot), FM (blue dashed), XY (pink dotted), and 
SF (green solid). In our simulation we use K = 20 for all 
phases and in the SF phase n\ — 112 — 1/2. For the XY 
phase 6a — 0b = 0.6-7T. 



V. CONCLUSIONS 

In summary we have generalized the model of a sin- 
gle component atomic lattice gas described by the Bose- 
Hubbard model coupled to near resonant optical cavities 
to the case of a two-component Bose-Hubbard model. 
We have shown conclusively through the probe cavity 
photon numbers and its spectra dependence on various 
system parameters that different quantum phases of the 
two-component Bose-Hubbard model can be easily dis- 
tinguished and confirmed. Our results shine new light on 
atomic lattice gases coupled to cavity QED systems. 
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